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ABSTRACT 

We have investigated the physical conditions under which accreting neutron stars in low-mass X-ray binaries 
can both produce and preserve sufficient quantities of carbon fuel to trigger superbursts. Our theoretical models 
span the plausible ranges of neutron star thermal conductivities, core neutrino emission mechanisms, and areal 
radii, as well as the CNO abundances in the accreted material. We find that neutron stars that accrete hydrogen- 
rich material with CNO mass fractions Zcno ^ Zcno.o will not exhibit superbursts under any circumstances. 
Neutron stars that accrete material with CNO mass fractions > 4Zcno,o will exhibit superbursts at accretion 
rates in the observed range. On this basis, we suggest that the mass donors of superburst systems must have 
enhanced CNO abundances. The accreted CNO acts only as a catalyst for hydrogen burning via the hot CNO 
cycle, and therefore it is the sum of the three elements' mass fractions, not the individual mass fractions 
themselves, that is important. Systems that exhibit superbursts are observed to differ from those that do not 
exhibit superbursts in the nature of their helium-triggered Type I X-ray bursts: the bursts have shorter durations 
and much greater a-values. Increasing the CNO abundance of the accreted material in our models reproduces 
both of these observations, thus once again suggesting enhanced CNO abundances in the mass donors. Many 
compact binary systems have been observed in which the abundances of the accreting material are distinctly 
non-solar. Though abundance analyses of the systems that exhibit superbursts currently do not exist, Bowen 
fluorescence blend profiles of 4U 1636-536 and Ser X-l suggest that the mass donor stars may indeed have 
non-solar CNO metallicities. More detailed abundance analyses of the accreting matter in systems that exhibit 
superbursts are needed to verify our assertion that the matter is rich in CNO elements. 

Subject headings: dense matter — nuclear reactions, nucleosynthesis, abundances — stars: neutron — X-rays: 
binaries — X-rays: bursts 



| Browil2004T:ICooper & NaravaiJ200i 120051 ICumming et alJ 
2005). Each of these models requires a carbon mass fraction 
> 0.1 at the base of the accreted layer to trigger a thermonu- 
clear instability. However, it is not understood why such a 
large amount of carbon should exist deep within the ocean. 
Previous theoretical calculati ons of both stea dy-state hydro- 
gen/helium burning (Schatzetal. 1999, 2003) and unstable 
helium ignition jSchatz et alJl200U 120031 iKoike et al]l2004t 
IWooslev et alJl2004t iFisker et alJ [2005) on the surface of an 
accreting neutron star yield far too little carbon to ignite a su- 
perburst. Consequently, researchers who model superbursts 
must set the value of the carbon mass fraction at the base 
of the accreted layer "by hand," with little physical motiva- 
tion. This is inadequate, since many superburst character- 
istics, such as the range of accretion rates in which super- 
bursts occur l|^ umming & Bildsten f200lUCooper & N aravan 
2005; Cumming et all 120051) and the amount of energy re- 
leased by neutrinos durin g a superburst JCumming & B ildsten 
2001; Strohmaver & Brown 2002; Cumming et al. 2005), are 
strong functions of the carbon mass fraction. Thus, no super- 
burst model can be fully self-consistent until the process by 
which the carbon fuel survives both stable and unstable burn- 
ing is understood. 

As of this writing, each system in which a superburst 
has been observed exhibits normal Type I X -ray bursts as 
weU JKuulkersI l200l Iu7t Zand et alJ l2004al) . The normal 
bursts these systems exhibit differ remarkably from the nor- 
mal bursts of systems that accrete at similar rates and from 
which no superbursts have been observed. First, the average 
e-folding decay times are smaller than those of normal bursts 

Electronic address: rcooper@cfa.harvard.edu, bmukhopa@cfa.harvard.edu, dsteeghi@MaialWdieffli^ rftyray^h@^fer?aTvara?eW urstS ' Second, for each 



1. INTRODUCTION 

Superbursts are energetic thermonuclear flashes that oc- 
cur on the surfaces of accreting neutron stars in low mass 
X-ray binaries. They are thought to be caused by unsta- 
ble carbon burning deep within the accreted layer. Su- 
perbursts distinguish themselves from their hydrogen- and 
helium-burning Type I X-ray burst counterparts (which we re- 
fer to as "normal" bursts) by their ~ 1000 times larger burst 
energies a nd ~ 1000 times longer recurrence times (fo r re- 
views, see Cumming 2004; Strohmaver & Bildsten 2005). As 
of this writing, nine superbursts have been observed in seven 
sources. Each of these superbursts has an integrated pho- 
ton flux of k 10 42 ergs and occurs in a system with an ac- 
cretion rate M » 0.1-0.2 5Mp,dd, where M^di denotes the Ed- 
dington limit (Kuulkers 2004). Four superburst candidates 
have been observed in GX 17+2, which has an accretion 
rate M sa QMIvm din't Zand et alJl2004albl) Very recently, 
iRemillard et al.1 ( 2005a) and Kuulkers (2005) have observed 
two likely superbursts in the systems 4U 1608-522 and 4U 
06 1 4+09 1 . The paucity of observational data makes the recur- 
rence times of superbursts difficult to determine, though three 
superbursts have b een observed wi t hin 4.7 years from the sys- 
tem 4U 1636-536 JWiinandsl2001t IStrohmav er & Markwardl 
2002; Kuulkers et all2004l) . 

Theoretical studies of superbursts have been rather success- 
ful at reproducing the general observational characteristics, 
such as the energetics, recurrence times, and absence of su- 
perbursts in systems with ac cr etion rates below sa . IMrm 
(ICumming & Bildstenl 1200 It IStrohmaver & Brownl 12001 
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system that has exhibited a superburst, a > 1000, where a 
is defined as the ratio of the energy released between normal 
bursts to the energy released during a normal burst. This is 
significantly greater than the a-values of most systems with 
similar a ccretion rates and in which no su perbursts h ave bee n 
detected (Ku ulkers et alJ2002btlin'tZand et alJ2003[ l2004b). 
For systems with high accretion rates and in which the ac- 
creted matter is predominantly hydrogen, iNaravan & Hevll 
(2003) found that normal bursts occur in a unique regime 
that they refer to as "delayed mixed bursts." In these sys- 
tems, a large fraction of the hydrogen and helium fuel burns 
stably before the full thermal instability is trigg ered. This 
stable b urni ng explains the observat ions of[£an Paradii set alJ 
( 1988aj and lCornelisse et all ( 120031) who found that, for sys- 
tems that accrete at a high rate and exhibit such Type I X- 
ray bursts, the quan t ity a r ises dramatically to values > 1000. 
Cooper & Naravan ( 2005) hypothesized that delayed mixed 
bursts may be the source of the substantial amount of carbon 
fuel needed to trigger a superburst in systems for which the 
accreted material is predominantl y hydrogen. Howev er, a- 
values of the delayed bursts of Narava n & Hevll (|2003) reach 
~ 1000 only if M ~ 0.3M Edd . For M ~ 0.1M Edd , roughly 
the lowest accretion rate at which delayed mixed bursts oc- 
cur, they find a < 100. The low a-values for M ~ 0.1M E dd 
are consistent with the results of previous and subsequent un- 
stable helium ignition studies, all of which yield a negligible 
amount of carbon. As of yet, no model has reproduced the 
observed characteristics of normal bursts from systems that 
exhibit superbursts. 

Any successful theoretical model of the low mass X-ray 
binaries considered in this paper must reproduce the follow- 
ing four phenomena that are unique to these systems. First, 
it must reproduce superbursts with energies and recurrence 
times that are consistent with observations. Second, nuclear 
burning in the oceans must produce sufficient amounts of 
carbon, and the carbon must survive until a superburst ig- 
nites. Third, delayed normal bursts with a > 1000 must occur. 
Fourth, the normal bursts that occur must have unusually short 
durations. Previous authors have made significant progress on 
the first phenomenon, but none has successfully addressed the 
latter three. In this investigation, we construct models that at- 
tempt to explain all four phenomena. We begin in |2] with 
a description of our theoretical model. In ^3] we discuss the 
physical conditions under which sufficient amounts of carbon 
may survive long enough to ultimately fuel a superburst. We 
discuss the results of our model in |4] and we compare our 
results to those of previous theoretical studies in |5] In 5j6] 
we compare our results with observations. Our model pre- 
dicts that the matter accreted onto the surface of the neutron 
star must be rich in CNO in order for superbursts to occur. 
We discuss the observational evidence that such an overabun- 
dance exists in and we conclude with a summary in |8] 

2. THE MODEL 

In this investigation, we assume that matter accretes spher- 
ically onto a neutron star of gravitational mass M and areal 
radius R at a rate M, where M is the rest mass accreted per 
unit time as measured by an observer at infinity. We set the 
hydrogen and heavy element composition of the accreted mat- 
ter to be that of the Sun, such that at the stellar surface the 
hydrogen mass fraction X out = 0.7 and the heavy element frac- 
tion Z out = 0.004, where Z refers to all metals other than CNO. 
We set the helium mass fraction Y out = 1 -X ut _ Z out -ZcNO,oub 
where ZcNO.out is the mass fraction of CNO elements accreted 



from the companion star. We treat Zcno.oui as a free parame- 
ter. 

2.1. Thermal and Hydrostatic Structure of the Crust 

To calculate the equilibrium configura tion of the neutron 
star cr ust, we use the theoretical model of Coope r*& Naravanl 
(2005) with two modifications. To calculate the inner tem- 
perature boundary condition at the crust-core interface, the 
authors assume that the core emits neutrinos via either modi- 
fied Urea (mUrca) reaction s or pionic reactions. They us e the 
neutrino luminosities from Shapiro & Teukolskv ( 1982b and 
thereby determine the core temperature. The core tempera- 
tures resulting from these neutrino luminosities are ~ 3 x 10 8 
K and ~ 2 x 10 7 K, respectively. However, if baryons ex- 
ist i n the core, they ar e probably superfluid at high densi- 
ties ( iBavm et alJl!969l) . Baryon superfluidity greatly sup- 
presses neutrino emission via mUrca reactions and there- 
fore raises the core tempe rature llYakovlev et aljfl99^.l200ll 
Yakovlev & Pethick 2004). Thus, when we model the effects 
of baryon superfluidity, we set the neutrino luminosity to be 



L i , = (2.0x 10 J7 ergss _1 ) 



Vio 9 k; 



(i) 



where m is the interior gravitational mass and T is the proper 
temperature, both evaluated at the crust-core interface. The 
core temperatures resulting from this neutrino luminosity are 
~6x 10 K. Note that we set the coefficient of L v not to corre- 
spond to any specific emissivity model, but rather to produce 
c ore te mperatu res in the desi red range. 

Cooper & Naravan (2005) define the free parameter Cf as 
the fraction of hydrogen and helium that ultimately burns to 
carbon. This parameter enables the authors to set the value 
of ZcNO,base, the mass fraction of carbon at the base of the 
accreted layer, such that 



ZcNO .base - ZcNO.out + Cf (X out + Tout) ■ 



(2) 



The value of Q should affect not only the carbon yield, but 
also the amount of energy generated within the crust, be- 
cause the total energy released per unit mass when helium 
is burned ultimately to iron peak elements is greater than the 
energy released when helium is burned only to carbon. To 
account for the additional energy generation, we include Cf 
in the energy conservatio n equation. Thus, equation (6) of 
Cooper & Naravan (2005) becomes 



d_ 

d£ V R 2 



1 + -Hl-C f ) 



(3) 



2.2. Nucleosynthesis in the Accreted Layer 

The thermal and hydrostatic structure model described 
above is quite sophisticated, approaching the state-of-the-art 
of such one-dimensional models. Its nuclear reaction net- 
work, however, is quite simplistic, since it includes reaction 
rates only for hydrogen, helium, and carbon burning. The hot 
nuclear flow through stable thermonuclear burning on the sur- 
face of an accreting neutron star almost certainly produces a 
wide assortment of different isotopes through hydrogen and 
triple-a reactions, rp- and ap-processes, and a-captures. To 
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calculate the detailed nuclear flow due to stable thermonu- 
clear burning, we consider a large number of isotopes on ei- 
ther side of the stability line. Our nuclear reaction network 
contains 255 nuclear species up to 72 Ge and all of the possible 
reactions between the various isotopes. Stable burning at the 
sub-Eddington accretion rates we consider in this investiga- 
tion produces negligible abundances of isotopes with atomic 
weights larger than 72 (e.g. lSchatz et all 199% . so the size of 
our network is adequate for our purposes. Lists of the ma- 
jor nuclear reactions t hat may take place in the flo w are given 
in previous work (e.g. Clavto Ml983llLang|ll999|). In our nu- 
cleosy nthesis code, we use the reaction rates of [Fo wler e^ flll 
(1975) but include the updated rate s of I Wagoner et alJ 
(119671). iFulleret all (119801 Il982albl). iThielemannl < 19801) . 
Walla ce & Wooslevl (119811) . and lHarris et all (11983k Each 
of our reaction rates has been updated and is current up 
to the year 1996 or later (F.-K. Thielemann 1996, private 
communication). The code has been successfully used by 
one of us in previous studi es of thermonuclear reactions 
in hot accretion discs jC hakrabarti & Mukhopadhvavj n"999t 
iMukhopadhvav & Ch akrabartil l2000l 120011) . Therefore, we 
are confident that this code is applicable to our present in- 
vestigation. 

The nuclear reaction network, which is a set of coupled dif- 
ferential equations, is linearized and evolved in time along 
the thermal and hydrostatic struc ture of the ocean derived 
from the model described in 32.11 T his well-proven m ethod 
is widely used in th e literature (e.g. Arnett & Trura rJll969t 
IWooslev et all 119731 KTukhopadh vav & Chakrabartil boOOl) 
Here, we briefly outline how we time-evolve the isotopic 
abundances. For simplicity, we consider only four iso- 
topes and only three reactions: l H(p,(3 + v)D, D(D,7) 4 He, and 
2 He(a,7) 12 C, although, in our actual calculations, we use 
255 isotopes and all of the possible reactions between the iso- 
topes as mentioned above. Neglecting any backward reac- 
tions, the corresponding rate equations can be expressed as 




/-A H 0^ 
A H -A d 
A D -A He 
V A He 0; 





/Xn\ 




Ad 


) 


Afle 







(4) 



where the various A's and A"'s are the reaction rates and mass 
fractions of the isotopes, respectively. The above equation can 
be written as 

dv 
dt 

If A is diagnolizable, we can solve this equation by finding 
the four eigenvalues A,- and four eigenvectors u,- of A. v(f ) is 
then a linear combination of the eigenvectors such that 



■■Ay. 



(5) 



v(f) = ^A,e A ''u,- 



(6) 



where one determines the constant coefficients A, by setting 
v(0) = Vq, the initial elemental abundance. Note that this 
method is correct only if the individual nuclear lifetimes are 
constant. In our actual computational code, we include 255 
isotopes, so A is a 255 x 255 matrix with 255 eigenvalues and 
255 eigenvectors, and thus i runs from 1 to 255. 

To evolve the composition of the entire layer in time, we 
perform this procedure at each timestep. The reaction rates 
are in general functions of the temperature, density, and com- 
position, and therefore they must be recalculated at each 
timestep. Normally, we set the timestep Af = AS/S, where 



AS is the thickness (mass per unit area) of a thin spherical 
shell and E is the mass accretion rate per unit area. However, 
we restrict the fractional change in the abundance of each 
species i to be less than 5, which we have set to 0.05 for our 
calculations. Specifically, for each species i such that its mass 
fraction^; > 10" 20 , we require Af < <5minK/(dX,/df)]- Note 
that the fractional change in the abundance of each species 
with a nontrivial mass fraction is always much less than 0.05. 
Nonetheless, we have performed calculations for values of S 
both somewhat greater than and much less than 0.05, and we 
have obtained the same results. Therefore, any error in our 
results due to our choice of stepsize is negligible compared 
to that due to the uncertainties in the reaction cross sections. 
We have verified that our nuclear reaction network conserves 
mass to within one part in 10 s . 

2.3. Carbon Abundance in the Accreted Layer 

To calculate the mass fraction of carbon that exists at the 
base of the accreted layer, we proceed as follows. First, we 
calculate the equilibri um configuration of the neutron star 
crust as described in 32.11 To perform this calculation, we 
choose an arbitrary value for the carbon mass fraction at the 
base of the layer. This value will be refined later. Second, 
we use this equilibrium configuration in our nucleosynthesis 
model to calculate the true carbon mass fraction at the base 
of the layer. In general, this value will be different than that 
used to determine the equilibrium configuration. Therefore, 
we repeat this two-step process, but we use the refined car- 
bon mass fraction derived from our nucleosynthesis code to 
calculate the equilibrium configuration. We iterate this pro- 
cess until successive values of the carbon mass fraction de- 
termined from our nucleosynthesis model agree to sufficient 
accuracy. This process converges rapidly, and the value to 
which it converges is completely insensitive to the value of 
our initial guess. As an example of a typical convergence, 
the carbon mass fraction after each iteration for a particular 
calculation, starting with an initial guess of 0.2, is 0.52709, 
0.56132, 0.56477, 0.56521, and 0.56531. After we have de- 
termined the carbon mass frac tion, we conduct the fu ll global 
linear stability analysis of Cooper & Naravan (2005) to deter- 
mine whether a superburst will occur, assuming that all of the 
carbon produced via stable burning has survived (see ^3}. 

The reader should note that we do not couple the nuclear 
energy generation rates of our reaction network directly into 
the energy conservation equation (0. In general, the contri- 
butions from hydrogen and helium burning dominate the total 
nuclear energy generation rate, and we calculate these con- 
tributions in equation (0 to high accuracy. Nevertheless, to 
ensure the internal consistency of our method, we have per- 
formed calculations in which the energy generation rates of 
our reaction network are coupled into equation (0, and we 
find that the differences in the final carbon mass fractions are 
insignificant. 

3. SURVIVAL OF THE CARBON FUEL 

Helium fusion via the triple-a reaction produces most of the 
carbon fuel that eventually triggers a superburst. Uncertain- 
ties in the carbon abundance deep within the accreted layer lie 
not in carbon's production, but in its survival. To ultimately 
become fuel for a superburst, carbon must survive both the 
stable burning that produces it and the unstable burning dur- 
ing normal Type I X-ray bursts that potentially consumes it. 

3.1. Survival During Stable Burning 
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The rapid proton (rp) process of Wall ace & Wooslev! ( I 19811) 
is the primary culprit in the destruction of 12 C. If hydro- 
gen is present, hydr ogen will burn via the hot CNO cy- 
cle (iHovle & Fowlerlfl 965ft . During the hot CNO cycle, es- 
sentially all of the CNO elements will be processed into 
14 and ls O. For the sub-Eddington accretion rates con- 
sidered in this paper, the dominant breakout reaction from 
the hot CNO cycle into the rp-process is ls O(a, 7) 19 Ne 
(Walla ce & Wooslevl U98TI iSchatz et alJ U999t IFisker et al.1 
120041) . It is possible for 19 Ne to return to the hot CNO cycle 
by the series of reactions I9 Ne(/3 + z/) I9 F(/;>, a) ie O(p, 7) 17 F(p, 
7) 18 Ne(/3 + i/) 18 F(p, a) 15 0, but if 19 Ne captures a proton by the 
reacti on 19 Ne(p, 7) 2Q Na, the ion can never return to the CNO 
cycle (Wallace & Wooslev 1981). Therefore, a 12 C ion that is 
either accreted from the companion star or produced from the 
triple-a reaction will not survive if it is both converted to ls O 
in the hot CNO cycle and removed from the hot CNO cycle 
via breakout reactions. 

For the relatively high accretion luminosities considered in 
this investigation, the rate at which hydrogen burns is set by 
the /3-decay timescales of 14 and 15 0, and therefore the rate 
is both temperature- and density-independent. The hydro- 
gen nuclear energy generation rate for the hot CNO c ycle is 
e H = 6 x 10 15 Z CN o ergsg -1 s _1 ( IHovle & Fowlerlll965k The 
column depth (mass per unit area) at which hydrogen burns 
out, Eh, is thus 

Sh «E^h« ( Sx1100s)^^, (7) 

CH ZcNO.out 

where £ is the mass accretion rate per unit area and is 
the total nuclear energy released per unit mass of hydrogen 
burned. In general, the column depth at which hydrogen burns 
increases with increasing accretion rate. 

In contrast, the rate at which helium burns is very sensi- 
tive to both temperature and density. The temperatures of the 
envelope and ocean of an accreting neutron star typically in- 
crease as the accretion rate increases. Consequently, for a pure 
column of helium, the column depth at which helium burns 
decreases with increasing accretion rate. This statement re- 
garding the depth at which helium burns is not strictly true 
for a hydrogen/helium mixture because the hydrogen burning 
will steadily add helium to the matter, but the general idea is 
still valid. 

If a large amount of carbon is to survive episodes of sta- 
ble burning to eventually trigger a superburst, most of the 
accreted hydrogen must burn via the hot CNO cycle before 
helium starts to burn. If this occurs, then the carbon produced 
via helium burning will not be processed into 15 and leave 
the CNO cycle through rp-process breakout reactions. Since 
the depth at which hydrogen burns increases with increasing 
accretion rate, and the depth at which helium burns decreases 
with increasing accretion rate, one would expect that more 
carbon wil l survive episodes of stab le burning at lower accre- 
tion rates dSchatz et al.lll999l l200l). Figure [T] illustrates this 
trend. 

3.2. Survival During Unstable Burning 

Any carbon that survives stable burning must also with- 
stand unstable burning in normal Type I bursts if the car- 
bon is to ultimately become fuel for a superburst. Pre- 
vious time-dependent models of normal prompt bursts, for 
which a < 100, imply that essentially no carbon will survive 
unstable hydrogen and helium burning dSchatz etaTll200lL 



2003| iKoike et al.1 120041 iWoosleveTaTl l2004t IFisker et alJ 
20051) . iCooper & Naravanl ( 120051) speculate that delayed 
bursts, which are bursts that are triggered after a long period 
of stable burning has occurred and thus result in extremely 
high a-values, will leave most of the thick layer of carbon 
that exists below the ignition region intact. Presumably, the 
increase in temperature due to hydrogen and helium burn- 
ing in the delayed burst is insufficient to ignite most of the 
carbon that exists deep in the ocean. Although we are cur- 
rently unable to either confirm or refute this assertion, we 
note that observations do imply that a sufficie nt amount of car - 
bon does indeed survi ve delayed bursts ( in't Zand et al. 2003). 
iWooslev et alJ ( 120041) show that the ashes of a given normal 
burst are reprocessed in subsequent bursts. This "composi- 
tional inertia" may destroy some fracti on of the carbon that 
a given burst leaves intact. However, Woo slev et alJ (120041) 
find that the effects of compositional inertia are diminished 
if the CNO abundance of the accreted material is high or if 
the burst recurrence time is long, both of which apply for the 
models we will suggest later in |4] Unfortunately, no detailed 
time-dependent studies of the delayed bursts observed in these 
systems have been carried out as of this writing. Therefore, 
our speculation that a large fraction of the carbon produced 
via stable burning will survive delayed bursts must be investi- 
gated further. For the purposes of this work, we assume that, 
if carbon is produced via stable burning, it will be destroyed 
if a normal prompt burst ignites (for which a < 100), and it 
will survive if either a normal delayed burst ignites (for which 
a 3> 100) or the system is stable to normal bursts. 

4. RESULTS 

We have constructed a total of 144 different models to deter- 
mine the conditions under which accreting neutron stars both 
produce and preserve sufficient amounts of carbon to trigger 
superbursts. We choose four different stellar areal radii, three 
different core neutrino emission mechanisms, and three differ- 
ent conductive opacity prescriptions, which likely bracket the 
true radii, core neutrino emissivities, and conductive opacities 
of neutron stars found in nature. See Table 1 for a list of these 
parameters. The column "Core v Emissivity" describes the 
neutrino emission mechanism in the core, where "hot mUrca" 
refers to a stellar core that cools via mUrca reactions sup- 
pressed by baryon superfluidity (see SI2.1> . "mUrca"' refers 
to a core that cools via mUrca reactions, and "Pion" refers to 
a core that cools via pionic reactions. The column "Q" de- 
scribes the conductive opacity of the crust. For this column, 
"5.2" and "100" refer to neutron stars with inner crusts that 
have formed orde red crystal lattices and which have im purity 
parameter values (Itoh & Kohvama 1993; Brown 2000) of 5.2 
and 100, respectively, and "disordered" re fers to a neutr on star 
with a completely disordered crust. See Brown (2004|) and 
Cooper & Naravan (2005) for details regarding these param- 
eters and their effects upon superbursts characteristics. Addi- 
tionally, we choose four different values of ZcNO.out, the mass 
fraction of CNO elements accreted from the companion star. 

A superburst will occur if a large amount of carbon deep 
within the ocean undergoes unstable thermonuclear fusion. In 
this investigation, we say that a superburst occurs if t he ca rbon 
produced by stable helium burning survives (see 33. Il l and 
the resulting carbon-rich accreted column is unstable accord- 
ing to the general-re lativistic global linear stability analysis of 
Cooper & Naravan (2005). For definiteness, we use the crite- 
rion that the carbon yield derived from our nucleosynthesis 
model will not survive if a normal burst occurs with a < 500. 
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FIG. 1.— 'H, 4 He, and 12 C mass fractions in the accreted column due to stable burning for two neutron stars with different accretion rates. M = 0. lMndd f° r 
Figure (a) and M = 0.3MEdd f° r Figure (b). Any 12 C accreted from the companion star is depleted when 1 H begins to burn because 12 C is processed into 14 and 
ls O during the hot CNO cycle. For M = O.lMEdd. 'H burns out before most of the 4 He has stalled to burn, so much of the 12 C produced via triple-a reactions 
survives. For M = 0.3MEdd, l H and 4 He burn simultaneously, so the 12 C produced via triple-a reactions is lost due to CNO cycle breakout reactions. 



This cutoff value is close to that of the superbursting system 
with the smallest a-value, 4U 1636-536, for which a ~ 440 
fcn't Zand et al.ll2004bl) . The precise value of this cutoff is 
unimportant for the final results. 

4.1. Results for Z C no,oui = Zcno,© 

For the relatively high accretion rates at which superbursts 
occur, the envelope and ocean of the neutron star are thermally 
insulated from the inner crust and core, and so the thermal 
profile of the ocean depends primarily on the accretion rate 
( Cooper & Naravan 2005). Consequently, the carbon yield 
resulting from stable hydrogen and helium burning is rather 
insensitive to both the conductive opacity of the crust and 
neutrino emission mechanism of the core, but it is quite sen- 
sitive to the accretion rate. We plot the 12 C mass fraction as a 
function of accretion rate for nine neutron stars with different 
conductive opacities and core neutrino emissivities in Figure 
|5Ji. The parameter Z acc = M/Me<h is the accretion rate normal- 
ized to the Eddington limit, where Msad = 4-irGM( l+z) /czK es , 
z= (1 — 2GM / Rc 2 )~ l I 1 — 1 is the gravitational redshift, and 
K es = 0.4 cm 2 g _1 . At lower accretion rates, most of the hy- 
drogen burns before the helium ignites, so the carbon pro- 
duced via stable helium burning survives. As the accretion 
rate increases, the hydrogen and helium ignition regions over- 
lap, so the carbon produced via stable helium burning will 
be processed into oxygen during the hot CNO cycle, and the 
oxygen will be rem oved from the hot CNO cycle via breakout 
reactions (see 33. 1> . Therefore, the carbon yield substantially 
decreases with increasing accretion rate. 

Unlike the conductive opacity and core neutrino emissiv- 
ity, the stellar radius significantly affects the carbon yield for 
a given / acc , as illustrated in Figure There are two rea- 
sons for this. First, the accretion rate per unit area E is lower 
for a larger radius, so Eh, the column depth at which hydro- 
gen burns, decreases with increasing radius. Second, the den- 
sity and temperature at a given column depth E are generally 
lower for a larger radius, so the column depth at which he- 
lium burns increases with increasing radius. The combination 
of these two effects implies that helium is more likely to burn 
in a hydrogen-deficient medium if the stellar radius is large. 
Therefore, stable burning on neutron stars with larger radii 
will yield more carbon fuel. 

The derived carbon yield for a given Z acc shown in Figure 



|2] is meaningless, however, if prompt normal bursts occur at 
that Z acc . We plot the a-values of normal bursts as a function 
of accretion rate for nine neutron stars with different conduc- 
tive opacities and core neutrino emissivities in Figure [3^ and 
with four different stellar areal radii in Figure [SJ). To per- 
form these normal burst calculatio ns, we use the model de- 
scribed by Remillard et al.l (12005b) tailored for persistent ac- 
cretors, which is a slightly modified version of the model of 
Cooper & Naravan (2005). Like the carbon yield due to sta- 
ble burning, a depends weakly on the thermal structure of the 
crust and core, but it is quite sensitive to the stellar radius. 
For the R = 10.4 km models shown in Figure^, no normal 
bursts occur for Z acc > 0.3. Delayed mixed bursts occur for 
0.2 < / acc < 0.3, prompt bursts occur for 0.05 < Z acc < 0.2, 
and delayed helium bursts occur for / acc < 0.0 5 (for details on 
the various bursting regimes, see|Naravan & Hevll2003l) . 

We list the range of accretion rates in which superbursts 
occur for each model in Table 1. No superbursts occur 
at any accretion rate for any of the 36 models for which 
ZcNO.out = Zcno,©, where Z CN o,o = 0.016. For the high ac- 
cretion rates at which either a> 100 or normal bursts do not 
occur, stable burning pro duces a negligible am ount of carbon, 
in agreement with lSchatz et all (119991 . 12003). For / acc ~ 0.1, 
prompt bursts destro y any carbon produ ced via stable burning 
(ISchatz et alJ 120011 120031 iKoike et all 120041 IWooslev et all 
2004: iFisker et al.ll2005l) . At still lower Z acc , the carbon may 
survive the delayed bursts, but the carbon that remains will 
burn stably ( Cumming & Bildsten 2001; Cooper & Naravan 
120051 ICummin get alJ2005l) and no superburst will occur. In- 
creasing the radius does not raise the probability that a super- 
burst will occur. Though increasing the radius will increase 
the carbon yield due to stable burning, it will also increase 
the limiting accretion rate above which delayed mixed bursts 
occur. These results hold for any ZcNO.out < Zcno,© as well, 
for decreasing ZcNO.out will both lower the carbon yield and 
extend the range of accretion rates over which prompt normal 
bursts occur. 

4.2. Results for Z C no,om > Z C no,q 
As mentioned in 33.11 carbon will survive to ultimately 
trigger a superburst only if the stable helium burning that 
produces it takes place in a hydrogen-deficient environment. 
Thus, the column depth at which the accreted hydrogen burns, 
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FIG. 2. — I2 C mass fraction at the the base of the accreted column due to stable burning as a function of accretion rate. The parameter Z acc is the accretion rate 
normalized to the Eddington limit. Figure (a) shows the 12 C mass fractions deep within the accreted layer of 10.4 km neutron stars with ZcNO.out = Zcno,0 for 
each of the nine thermal structure models. Figure (b) shows the 12 C mass fractions of neutron stars with ZcNO.out = Zcno.0> Q = 5.2 below the accreted layer, 
and with cores that emit neutrinos via mUrca reactions for four different radii. 




Eh, must be less than the column depth at which helium burns 
to provide enough fuel for a superburst. According to equa- 
tion Q, Eh is inversely proportional to ZcNO.out, so if the 
accreted matter is rich in CNO elements, hydrogen is more 
likely to burn out before helium ignites. Therefore, more of 
the carbon produced via stable helium burning should survive. 
Figure[4]illustrates this effect of the CNO mass fraction of the 
accreted matter on the carbon yield. We emphasize that the in- 
dividual mass fractions of the accreted carbon, nitrogen, and 
oxygen are unimportant. Since the accreted CNO ions are 
simply catalysts for hydrogen burning, it is only the sum of 
the three indi vidual mass fractions that matters. 

As in 34.11 the derived carbon yield for a given / acc is mean- 
ingless if prompt normal bursts occur at that Z acc . We plot 
the a-values of normal bursts for four neutron stars with dif- 
ferent CNO abundances of the accreted matter in Figure [5] 
Evidently, increasing ZcNO.out affects normal bursts in several 
ways. First, it lowers the critical accretion rate above which 
normal bursts do not occur. Hydrogen burning is always j3- 
limited at the accretion rates considered in this article, so only 
unstable helium burning triggers normal bursts. The helium 
nuclear energy generation rate eee monotonically increases 
with temperature, but the temperature sensitivity decreases 
as temperature increase s. In other words, d 2 e^ e /dT 2 < 
(Fushiki & Lambl U987l) . Therefore, helium burns stably at 
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FIG. 4. — l2 C mass fraction at the the base of the accreted column due 
to stable burning as a function of accretion rate for four different values of 
ZcNO.out. the mass fraction of CNO elements in the accreted matter. R = 10.4 
km, Q = 5.2 below the accreted layer, and the stellar core is assumed to emit 
neutrinos via mUrca reactions. 



sufficiently high temperatures (Bildsten 1998). Since the tem- 
perature of the ocean increases with accretion rate, there exists 
a critical accretion rate above which normal bursts do not oc- 
cur. The increased hydrogen nuclear energy generation rate 
due to the CNO enhancement raises the temperature of the 
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TABLE 1 
Results 



M/Mrjdd Range for Superbursts for Zcno.oui = 



ft 


Lore v 




^CNO,0 2ZcNO,Q 


J^CNO,© 


4A:no,0 


[Km) 


ninissiviiy 




8.2 


hot mUrca 


5.2 








8.2 


hot mUrca 


100 




0.04-0.06 


0.04-0.08 


8.2 


hot mUrca 


disordered 




0.03-0.07 


0.03-0.10 


8.2 


mUrca 


5.2 




0.05-0.07 


0.03-0.09 


8.2 


mUrca 


100 




0.03-0.07 


0.03-0.11 


8.2 


mUrca 


disordered 




0.03-0.07 


0.03-0.13 


8.2 


Pion 


5.2 




0.03-0.07 


0.03-0.06 


8.2 


Pion 


100 




0.03-0.07 


0.03-0.11 


8.2 


Pion 


disordered 




0.03-0.07 


0.03-0.14 


10.4 


hot mUrca 


5.2 


0.04-0.05 


0.06-0.09 


0.07-0.13 


10.4 


hot mUrca 


100 


0.03-0.05 


0.06-0.09 


0.07-0.13 


10.4 


hot mUrca 


disordered 


0.03-0.05 


0.06-0.09 


0.07-0.13 


10.4 


mUrca 


5.2 


0.03-0.05 


0.06-0.09 


0.07-0.13 


10.4 


mUrca 


100 


0.03-0.05 


0.06-0.09 


0.07-0.13 


10.4 


mUrca 


disordered 


0.03-0.05 


0.06-0.09 


0.07-0.13 


10.4 


Pion 


5.2 


0.03-0.05 


0.06-0.09 


0.07-0.13 


10.4 


Pion 


100 


0.03-0.05 


0.06-0.09 


0.07-0.13 


10.4 


Pion 


disordered 


0.03-0.05 


0.06-0.09 


0.07-0.13 


13.1 


hot mUrca 


5.2 






0.10-0.18, 0.22 


13.1 


hot mUrca 


100 






0.10-0.18, 0.22 


13.1 


hot mUrca 


disordered 






0.10-0.18, 0.22 


13.1 


mUrca 


5.2 






0.10-0.18, 0.22-0.25 


13.1 


mUrca 


100 






0.10-0.18, 0.22 


13.1 


mUrca 


disordered 






0.10-0.18, 0.22 


13.1 


Pion 


5.2 






0.10-0.18, 0.22-0.25 


13.1 


Pion 


100 






0.10-0.18, 0.22-0.25 


13.1 


Pion 


disordered 






0.10-0.18, 0.22-0.25 


16.5 


hot mUrca 


5.2 






0.04-0.07, 0.14-0.32 


16.5 


hot mUrca 


100 






0.04-0.07, 0.14-0.32 


16.5 


hot mUrca 


disordered 






0.04-0.07, 0.14-0.32 


16.5 


mUrca 


5.2 






0.04-0.07, 0.14-0.32 


16.5 


mUrca 


100 






0.04-0.07, 0.14-0.32 


16.5 


mUrca 


disordered 






0.04-0.07, 0.14-0.32 


16.5 


Pion 


5.2 






0.04-0.07, 0.14-0.32 


16.5 


Pion 


100 






0.04-0.07, 0.14-0.32 


16.5 


Pion 


disordered 






0.04-0.07, 0.14-0.32 



"For some models, superbursts occur over two distinct ranges of M/MecM separated by a stable 
zone. Both M/MecM ranges are given in these cases. 



ocean, so helium is more likely to burn in a stable fashion, re- 
sulting in a lower critical accretion rate. Second, the a-values 
of normal bursts are generally higher. Just below the criti- 
cal accretion rate, the normal bursts that occur are "mixed" 
bursts, meaning that the bursts consume substantial amounts 
of both hydrogen and helium. Since Sh is lower, a larger 
fraction of the hydrogen burns before helium ignites unstably, 
which means that less hydrogen is consumed during the burst. 
This results in a higher a. Third, the region of delayed helium 
bursts is curtailed, and a range of accretion rates exists below 
this region in which normal bursts do not occur. Again, this 
is due to the increased temperature of the ocean due to the 
intensified hydrogen energy generation rate. 

We list the ranges of accretion rates at which superbursts 
are triggered for Z CN o.out = 2Z CN o,©, 3Z CN o,o, and 4Z CN o,© 
in Table 1 . We find that an enhancement of the CNO abun- 
dance in the accreted matter will allow for superbursts to oc- 
cur in certain ranges of accretion rates. The quantity ZcNO.out 
is unique in that it is the only free parameter that will both 
increase the carbon yield due to stable burning at a given ac- 
cretion rate and decrease the range of accretion rates at which 
prompt normal bursts occur. Therefore, we suggest that an 
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?CNO = 0064 
J Z CN0 = 0.048 
Z CN0 = 0.032 



3NO 



= 0.016 
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FIG. 5. — a-values of normal bursts as a function of accretion rate for four 
different values of ZcNO.out- R = 10.4 km, Q = 5.2 below the accreted layer, 
and the stellar core is assumed to emit neutrinos via mUrca reactions. 



enhancement of the CNO abundance in the accreted matter is 
a prerequisite for superbursts to occur in systems that accrete 
predominantly hydrogen. 
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5. COMPARISON WITH PREVIOUS THEORETICAL STUDIES 

As noted in 34. II our conclusion that neutron stars that 
accrete material with ZcNO.out ;$ Zcno.o W1 U not exhibit su- 
perbursts in consistent with the results of previous theoreti- 
cal investigations. However, a comparison between Figure^ 
of this work and Figure 2 of lSchatzetail (^003) illustrates 
that a discrepancy exists between the mass fraction of car- 
bon that survives stable burning as a func tion of accretion rate 
that we calculate and the same quantity [Schatz etafl (120031) 
derive. While both groups find that the carbon mass frac- 
tion » 0.3 at l arc = 0.1 (t he lowest accretion rate they con- 
sider), Schatz et al. (2003) find that the mass fraction w 0.08 
at Z acc = 0.3, while our calculations imply that a negligible 
amount of carbon will survive at this accretion rate. Gener- 
ally speaking, it appears that the resulting carbon mass frac- 
tion is more sensitiv e to the accretion r ate in our model than 
it is in the model of Schatz et al. (2003). In particular, Figure 
[Q) illustrates that the hydrogen mass fraction plateaus around 
S 10 8 gem" 2 in our model, whe reas no such plateau exists 
in the model of Schatz et al. (2003). In this section we discuss 
our efforts to resolve this discrepancy. 

One difference between the two models is that ours is 
general-relativistic, while theirs is Newtonian. Consequently, 
the gravitational acceleration in our model is greater. Al- 
though this will reduce the carbon yield somewhat, the ef- 
fect is too small to account for such a large disparity. Fur- 
thermore, we would expect that the larger gravitational accel- 
eration would reduce the carbon yield at all accretion rates, 
which means that the results of the two models should dis- 
agree at all accretion rates. Thus, we deduce that the inclu- 
sion of general relativity in our model is not the source of the 
discrepancy. 

Another difference between the two models is the treat- 
ment of the inner boundary condition. While we set the inner 
bound ary condition at the crust-core interface, Scha tz et"aT1 
(2003) set it at the base of the burning layer. Therefore, we 
have conducted experimental calculations in which we set our 
boundary condition near the base of the burning layer. We are 
unable to reproduce their results at / acc = 0.3 for any reason- 
able choice of boundary condition, and so we conclude that 
the treatment of the inner boundary condition is not the sole 
cause of the discrepancy. 

A potential deficiency in our model may be that we do not 
directly couple our reaction network into our energy conserva- 
tion equation. This could result in an inaccurate thermal pro- 
file which would affect our carbon yield. As discussed in 32.31 
we have carried out calculations, specifically at Z acc = 0.3, in 
which the nuclear energy generation rate of our reaction net- 
work is directly coupled into equation Q, and the final carbon 
yield differs negligibly from our old result. Different radiative 
opacity prescriptions would affect the thermal profile too, but 
we use the same prescription as Schat zetall ll2003l) . so the 
opacity is likely not an issue either. 

A plausible explanation for the discrepancy which we are 
not able to rule out is that some of the corresponding re- 
actions rates in the two networks differ significantly. The 
two networks are of similar vintage, however, so we pre- 
sume that none of the corresponding reaction rates disagrees 
considerably. Unfortunately, since we do not have access to 
the network of ISchatz et alJ {2003), we are unable to con- 
duct a detailed comparison. However, two particular reac- 
tions that would certainly affect the final carbon yield and 
the hydrogen mass fraction profile are the CNO breakout re- 
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FIG. 6. — 'H mass fraction profiles in the accreted column due to sta- 
ble burning at M = 0.3MecM for all possible combinations of the most re- 
cent 15 0(o,7) 19 Ne and 14 0(a, p)' 7 F reaction rates. The bottommost profile, 
where *H is depleted at E M 10 7 8 gem -2 , is the result of a calculation in 
which we artificially lowered both reaction rates by a factor of 10 and 100, 
respectively. 



actions 15 Q(a, V) 19 Ne and u O(a, p) 17 F ISchatz et alJI!999t 
Fiske r et al ] l200l . To investigate their effects, we repeated 
the calculation shown in Figure^) using all possible combi- 
nations of the three most recent 15 0(a, 7) 19 Ne rates and the 
two most recent 14 0(a, p) 17 F rates (J. L. Fisker 2005, pri- 
vate communication). Figure[6]shows the resulting hydrogen 
mass fraction profiles. The CNO breakout reactions extract 
14 and ls O from the CNO cycle, causing the hydrogen burn- 
ing to stall. Hydrogen thus survives to a greater depth, which 
reduces the carbon yield. Although the lower ls O(a, 7) 19 Ne 
and 14 0(a, p) 17 F reaction rates curtail the plateau somewhat, 
we find that none of the six possible reaction rate combina- 
tions eliminates the plateau entirely. Furthermore, all of the 
final carbon yields differ by less than a factor of 2. We are able 
to remove the hydrogen plateau only by artificially reducing 
the 15 0(a, 7) 19 Ne and 14 0(a, p) 17 F reaction rates by a factor 
of 10 and 100, respectively. When we do this, we find that the 
carbon mass fraction « 0.08 for Z acc = 0.3. 

In summary, we are unable to pinpoint the reason why the 
final carbon yield due to stable burning at relatively high ac- 
cretion rate s derived from ou r model differs from the yield 
derived by Sch atz et"afl ll2003l) . Fortunately, our general con- 
clusions remain unchanged regardless of which model is a 
better description of the nucle ar reactions tha t occur in nature, 
for even if the model of Scha tz et all (|2003) is "correct," we 
still find that the mass donors in systems that exhibit super- 
bursts must be evolved stars with enhanced CNO abundances 
in order to both generate sufficient amounts of carbon fuel and 
reproduce the observed characteristics of the normal bursts. 

6. COMPARISON WITH OBSERVATIONS 

The nine superbursts that have been observed as of this writ- 
ing, excluding those from the anomalous system GX 17+2, 
occurred in systems with accretion rates between 10% and 
25% of the Eddington limit. The normal bursts that these 
systems exhibit differ in two ways from the normal bursts 
in other systems with similar accretion rates. First, the nor- 
mal bursts are delayed bursts, with a > 1000 ( iKuulkers et alJ 
2002a; in't Zand et al1 l2003l l2004bl) . whereas most systems 
accreting at these rates exhibit prompt bursts, wit h a < 100 
( Ivan Paradiis et alJll988aL iNaravan & Hevll 120031) . Second, 
the average durations of normal bursts that occur in these 
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systems are shorter than those of normal bursts in system s 
that do not exhibit superbursts (lin't Zand et al.ll2003l l2004b). 
Any successful theoretical model of superbursts must be able 
to explain these facts. In this investigation, we concen- 
trate on systems in which the accreted matter is predomi- 
nantly hydrogen. This excludes the system 4U 1820-303, in 
which the accreted matter is probably d ominated by helium 
JStrohmaver & Brownl 120021 [Cummins! 120031) . It is likely 
that the long periods of stable helium burning that take place 
in the hydrogen-deficient ocean of this source produce large 
amounts of carbon fuel for superbursts. We note, however, 
that it is not understood theoretically why this system does not 
exhibit normal bursts during these periods, when the accre- 
tion r ate is near its maximum (Bilds tenll995l[l998tlCum ming 
12001 . 

We begin by discussing the range of accretion rates in 
which superbursts have been observed, / acc ks 0.1-0.25. From 
Table 1, we require that ZcNO.out ^ 4Zcno,o in order for a su- 
perburst to occur in this range. ZcNO.out is the most important 
parameter that determines whether a system in this range of 
accretion rates will exhibit a superburst. If the CNO mass 
fraction of the accreted matter is low, we find that there is no 
scenario in which a sufficient amount of carbon fuel will sur- 
vive to ultimately trigger a superburst. The second most im- 
portant parameter is the stellar areal radius. Although a neu- 
tron star with almost any plausible radius can exhibit a super- 
burst at Z acc w 0.1 given a sufficiently large value of ZcNO,out, 
neutron stars will not exhibit superbursts at Z acc ps 0.2 unless 
the radius R > 13 km. If all neutron stars are indeed rather 
large, then our model predicts that superbursts will not occur 
in systems that accrete predominantly hydrogen if Z acc < 0.1, 
in agreement with observations. Our model may also explain 
why superbursts are not observed in most low mass X-ray 
binaries with / acc > 0.3 (excluding the anomalous GX 17+2, 
which we will discuss below). 

Previous theoretical models of normal bursts all produce 
prompt bursts at / acc «0.1,witha< 100. The results of these 
models are inconsistent with the delayed bursts observed at 
Z acc w 0.1 from the systems 4U 1636-536, KS 1731-260, and 
4U 1254-690, in whic h a » 44 0, 780, and 4800, respectively 
llin't Zand et alJ 120031 | 2Q 04b). These models usually as- 
sume t hat ZcNO.out ^ Zcno.o- However, Tillett & MacDonald 
( 1992) showed that theoretical models of normal bursts from 
the superburster 4U 1636-536 are inconsistent with observa- 
tions unless a significant amount of CNO exists in the ocean. 
Figure [5] shows that an enhancement of the CNO abundance 
in the accreted matter produces delayed bursts at Z acc w 0.1 
that are consistent with observations. The systems Ser X-l, 
GX 3+1, and 4U 1735-444 h ave accretion rates / acc w 0.2, 
0.2, and 0.25 (Kuulkers 2004), and they exhib it normal bursts 
with a ss 580 0, 2100, and 4400, respectively (lin't Zand et al.1 
120031 l2004bl) . Figure |5] shows that delayed bursts occur at 
these accretion rates only if the CNO abundance in the ac- 
creted matter is high. Therefore, our theoretical models of 
normal bursts are consistent with the observed delayed bursts 
with a > 1000 only if the CNO abundance in the accreted 
matter in notably greater than solar. 

The average decay times of normal bursts in all systems that 
exhibit superbursts are lower than the average decay times 
of normal bursts in systems that do not exhibit superbursts 
llin't Zand et alJl2003l l2004bl) . We plot the effective normal 
burst duration fe+He as a function of Z acc in Figure [7] where 
?H+He equals the energy released via hydrogen and helium 
burning during a normal burst divided by the Eddington lu- 
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FIG. 7. — Effective normal burst duration as a function of accretion rate, 
calculated by dividing the total burst energy by the Eddington luminosity, for 
four different values of ZcNO.out- R = 10.4 km, Q = 5.2 below the accreted 
layer, and the stellar core is assumed to emit neutrinos via mUrca reactions. 
Note that ?H+He decreases with increasing ZcNO.out f° r 'acc 0.06. 



minosity (see Naravan & Hevl 2003). FigureQillustrates that 
systems with a large CNO abundance in the accreted matter 
will exhibit normal bursts with smaller durations. Again, our 
models of normal bursts are consistent with observations only 
if the CNO abundance in the accreted matter in high. 

Several of the systems have a-values that are even higher 
than those plotted in Figure [5] This could easily occur if 
the systems have accretion rates close to a critical accre- 
tion rate that separates s table burning and unstable burning 
Jvan Paradiis et alj fl979). The systems 4U 1636-536, KS 
1731-260, Ser X-l, GX 3+1, and 4U 1735-444 all undergo 
episodes of irregular bursting behavior, which implies that 
they are accreting near such a critical rate tLewin et all 19871: 
Muno et alJ2000t |Sztaino et a l.ll983tlden Hartog et alJ2003t 
van Paradii set alJl!988bl) . Previous theoretical models pre- 
dict that normal bursts will occur for all accretion rates be- 
low the critical r ate ab ove which helium burning is stable. 
iNaravan & Hevll (12003) find that the critical accretion rate is 
roughly 30% of the Eddington limit, which is consistent with 
observations (Cornelisse et al. 2003; Remillard et al. 2005b). 
Thus it is not surprising that Ser X-l, GX 3+1, and 4U 1735- 
444 experience this irregular behavior, for their accretion rates 
are close to this upper bound. What is surprising is that 4U 
1636-536 and KS 1731-260 experience irregular bursting be- 
havior too, since Z acc w 0. 1 for these two systems. This implies 
that the range of accretion rates at which normal bursts occur 
is not continuous, for there exists a range of accretion rates 
below Z acc ps 0.1 in which normal bursts do not occur. The 
bursting behavior of KS 1731-2 60, which is a tra nsient sys- 
tem, is particularly interesting. iMuno et all J2000I) find that 
short, photospheric radius expansion bursts occur at high ac- 
cretion rates, where / acc ~ 0.1, and long bursts with no evi- 
dence of radius expansion occur at low accretion rates, where 
Lcc ~ 0.01, but no bursts occur at intermediate accretion rates. 
All of this is nicely reproduced in Figures El and |71 fo r mod- 
els wit h CNO enhancement in the accreted matter. Shi h et alJ 
(2005) suggest that the accretion rate of 4U 1636-536 is cur- 
rently decreasing with time. Our models predict that in this 
system in the near future either a will rise or normal bursts 
will cease altogether. This prediction should be testable with 
further observations. 

The system GX 17+2, which has an accretion rate / acc w 0.8 
and exhibits normal bursts with a ~ 1000 as well as super- 
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bursts, is still a mystery to us. Though there is little di rect evi- 
dence that the accretion rate is so high, in ' t Zand et al.l lf2004a) 
conclude with reasonable certainty that the accretion luminos- 
ity is always greater than 60% of the Eddington limit. If this 
is indeed true, then clearly the observations of GX 17+2 are 
inconsistent with our results. Theoretical models predict that 
systems that accrete at this rate will show su perbursts given 
a sufficient amount of carbon l|C umming & Bildsten] 1200 It 
lBrownll2004t iCoooer & Naravanl 120051) . However, it is not 
understood why so much carbon should survive or why nor- 
mal bursts should occur at such a high accretion rate. We 
offer two possible explanations, neither of which is without 
issues. First, we find that a reduction of the radiative opac- 
ities by roughly an order of magnitude will produce delayed 
mixed bursts at l acc ps 0.8 with a ~ 1000 and substantially in- 
crease the carbon yield due to stable burning. GX 17+2 has a 
magnetic field that is much stronger than th ose of typical low 
mass X-ray binaries ( Kuulker s et alJ2002ah . and strong inter- 
nal radial magnetic fields lower the radiative opacities of the 
ocean (Meszaros 1992; Lamb 2000). However, if an internal 
magnetic field were solely respons ible for the opa city reduc- 
tion, it would have to b e ~ 10 13 G (Ivan R iper 1988), whereas 
Wiina nds et all JT996T) derive an upper limit of « 5 x 10 9 G 
for this system. Therefore, we think this explanation is highly 
unlikely. Second, we find that lowering the hydrogen mass 
fraction of the accreted matter, X out , by a factor of ~ 2 will 
produce delayed mixed bursts at l dcc » 0.8 with a ~ 1000 
and substantially increase the carbon yield due to stable burn- 
ing. This will lower the burst durations as well. This ex- 
planation is definitely plausible, for theoretical evolution- 
ary models of intermediate-mass X-ray binaries sugg est that 
the seconda ries will be hydrogen-poor ( Podsiadlows ki et all 
I27M 12001 . Though some of the normal bursts observed in 
GX 17+2 were quite short (^10 seconds), most of the ob- 
served normal bursts were rather long (^10 minutes), which 
implies that hydrog en was abundant in the ocean when the 
bursts occurred l|K ahn ffi: Grindlayll 19§4[ iTawara et alJll984l 
ISztaino et aflll986t iKuulkers et alJll997l l2002al) . This issue 
should be investigated further. 

7. OBSERVATIONAL EVIDENCE OF CNO ENHANCEMENT 

Evidence that the accreted material in compact binaries can 
have non-solar abundances is found in a variety of systems. 
For example, UV and X-ray spectroscopy of accreting white 
dwarfs has revealed anomalous N/C emission l ine ratios in the 
accret ed material in a number of system s (e.g.lGansicke et alJ 
120031 iBonnet-Bidaud & Mouchellbo'oH iRamsav et alJl2005l 
and references therein). Since the abundances of the accreting 
material directly reflect the abundances of the mass donor star, 
various evolutionary scenarios have been proposed to explain 
these significant deviations from solar abundances. Given 
that hydrogen-rich donor stars will produce CNO deep within 
their cores, CNO processing as the mass donor star burns hy- 
drogen along its main sequence will naturally produce CNO- 
enriched gas, with most of the C converted into 13 N during 
the CNO cycle due to the long /3-decay timescale of 13 N. If 
these products reach the surface, we may expect to see signif- 
ica nt non- solar abundances i n the accretion flow. For exam- 
ple, EhaiilinsiniEin J2Q02) find an expected (N/N Q ) ps 5-8 
if the mass donor star in the cataclysmic binary QZ Ser had 
undergone significant hydrogen burning before mass trans- 
fer onto the white dwarf had started. Alternatively, the mass 
donor may have lost (part of) its hydrogen envelope, expos- 
ing its CNO core and significantly boosting the CNO mass 



fraction of the accreting matter. Enhanced CNO abundances 
have a lso been report ed for the black-hole X-ray binary XTE 
J 1 1 1 8 (lHaswell et all2002l) . and lJimenez-Garate et all ( 2005 ) 
report 4 < (N/O)/(N/O) < 9 in the high-mass X-ray bi- 
nary Her X-l. Other possible, though perhaps less likely, 
avenues by which the mass donor stars can become CNO- 
enriched include carbon accretion onto the donor star during 
the asy mptotic giant branch phase of the neutron star pro- 
genitor (de Koo l & Greenll995l: ISteinhardt & Sasselo\ l20()5T) 
and ejecta capt ure from the su pernova that begot the neutron 
star (e.g. Israelian et al. 1999). Not many abundance anal- 
yses for X-ray binaries have been published, but it appears 
that the accretion of CNO-enriched material occurs in a num- 
ber of compact binaries harboring evolved mass donor stars. 
We n ote that that the soft X-r ay transient 4U 1608-522, in 
which Remillar d et alJ (|2p05a) observe d a likely superburs t, 
may contain an evolved secondary star ( Wachter et al. 2002). 

Unfortunately, abundance analyses for the systems that ex- 
hibit superbursts have not yet been reported. The large num- 
ber of strong C, N, and O spectral lines in the UV regime 
make it the best window for such a study, but these analyses 
are severely impeded by interstellar extinction at those wave- 
lengths. Strong H and He lines tend to dominate the opti- 
cal spectroscopy, and the fact that these systems persistently 
accrete at 0.1 < Z acc < 0.25 hampers the study of the pho- 
tospheric composition of the donor star since the accretion 
light dominates. For such systems, the Bowen fluorescence 
blend near A4640-4660 A has proven to be a useful in direct 
probe of the mass donor star (Steeghs & Casares 2002). This 
blend consists of several N III components that are part of the 
Bowen fluorescence process, as well as nearby C III lines. Ir- 
radiation of the mass donor star leads to sharp and resolvable 
emission components throughout this blend originating from 
its surface. This was first demonstrated in Sco X-l, where 
the N and C components from th e donor were detected a t 
similar strengths ( see Figure 1 in Steeghs & Casares 2002). 
iHvnes et alJ J2004) report optical spectroscopy of the super- 
burst source Ser X-l and remark that only the N III lines were 
detected in the Bowen emission blend, and no C III emission 

^^^^^^^^^^^^^^^ ^^^^^T. 

was evident. Casares et al. (2004) present a compilation of 
Bowen blend profiles from a number of X-ray binaries includ- 
ing the superburst sources V801 Ara (also known as 4U 1636- 
536) and V926 Sco (also known as 4U 1735-444). Again, in 
V801 Ara sharp components are detected from the N III com- 
ponents, but no sharp C III components are evident. For V926 
Sco, a spectral feature coinciding with one C III line is de- 
tected, but the second C III line is not detected making it less 
clear-cut for this source. 

We remark that, since the N III transitions are part of the 
Bowen fluorescence process, whereas the C III are not, this 
blend cannot be used for quantitative abundance analysis. 
However, the dominance of N over C in two of the three ob- 
served superburst systems is certainly suggestive that the ac- 
creting matter may indeed be CNO-enriched. A more detailed 
abundance analysis of these systems could test whether the 
composition of the accreting matter does in fact conform to 
the expectations of our model calculations. 

8. SUMMARY AND CONCLUSIONS 

By merging t he theoretical superburst model of 

Cooper & Naravan ( 200 ^ with the detailed nuclear re- 
action network of Mukhopadhvav & Chakrabarti (2000), we 
have investigated the physical scenarios in which accreting 
neutron stars can both produce and preserve sufficient 
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amounts of carbon fuel to trigger superbursts. We have con- 
structed a total of 144 different models that span the possible 
ranges of neutron star thermal conductivities, core neutrino 
emission mechanisms, and areal radii, as well as the CNO 
abundances in the accreted material. We find that neutron 
stars that accrete hydrogen-rich material with CNO mass 
fractions less than or equal to that of the Sun will not exhibit 
superbursts, regardless of their accretion rates, conductivities, 
core neutrino emissivities, or radii. Neutron stars that accrete 
material with CNO mass fractions > 4Zcno,o can exhibit 
superbursts at accretion rates in the observed range, but only 
if the stellar radii are sufficiently large. We remark that the 
accreted carbon, nitrogen, and oxygen serve only as catalysts 
for hydrogen burning, and therefore only the sum of their 
individual mass fractions is important, not the individual 
mass fractions themselves. 

Observationally, systems that exhibit superbursts differ 
from systems that do not exhibit superbursts not only in the 
amount of carbon that survives deep in the oceans of the neu- 
tron stars, but also in the nature of their normal bursts. Com- 
paratively, the normal bursts observed in systems that exhibit 
superbursts are much shorter in duration, and they have en- 
values that are much greater. Increasing the CNO abundance 
of the accreted material both decreases the durations of nor- 
mal bursts and increases their ct-values, in agreement with ob- 
servations. This is the first study in which the behavior of nor- 
mal Type I X-ray bursts has been incorporated into models of 
systems that exhibit superbursts. Observers have discovered 
many compact stellar X-ray sources in which the accreted ma- 
terial is significantly non-solar. However, abundance analyses 
of the accreted material in systems that exhibit superbursts 
currently do not exist. We suggest that the secondaries of the 
low-mass X-ray binary systems Ser X-l, 4U 1254-690, 4U 
1735-444, GX 3+1, KS 1731-260, and 4U 1636-536 are all 
rich in CNO elements. Although there is some indirect evi- 



dence in a few of these systems that the mass donor stars have 
indeed undergone CNO-processing, more observations are re- 
quired to either verify or refute this assertion. 

One issue we have not addressed is the disparity between 
the superburst energies and recurrence times that have been 
observed and the energies and recurrence times derived from 
theoretical models. Specifically, all models predict that the 
energies and recurrence times should be roughly an order of 
magnitude larger than those observed. Theory and observa- 
tions still disagree even if the thermal conductivit y of th e crust 
is low and the core temperature is high. Brown (2004) notes 
that observations of KS 1731-260 imply that the neutron star 
in this system in fact has a hig h thermal conductivity and a 
low core temperature dWiinands et alJ20"02i) . which illustrates 
the sev erity of this discrepancy. Furthermore. lCumming et all 
(2005) show that neutrino emission via Cooper pairing of su- 
perfluid neutrons in the inner crust limits the temperature of 
the ocean to values well below the w 6 x 10 8 K needed to 
trigger superbursts that match observations, and they suggest 
that an additional heating mechanism is required. Our mod- 
els do not provide for such a mechanism, for we find that any 
additional heating due to the increased CNO abundance of 
the accreted material has a negligible effect on superburst en- 
ergies and recurrence times. Superbursts may potentially be 
excellent probes into the interiors of neutron stars. Further 
progress in the physics of superbursts, combined with better 
observations, may ultimately lead to better understanding of 
not only neutron star interiors, but also fundamental physics. 

We would like to thank Jorge Casares, Paul Green, Josh 
Grindlay, Hendrik Schatz, Jeffery McClintock, and Charles 
Steinhardt for helpful discussions, Jacob Fisker for kindly 
providing updated reaction rates, and the referee for useful 
comments and suggestions. D. S. acknowledges a Smithso- 
nian Astrophysical Observatory Clay Fellowship. This work 
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